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Entanglement measures constructed from two positive, but not completely positive maps 
on density operators are used as constraints in placing bounds on the entanglement of for- 
mation, the tangle, and the concurrence of 4 x N mixed states. The maps are the partial 
transpose map and the $-map introduced by Breuer [H.-P. Breuer, Phys. Rev. Lett. 97, 
080501 (2006)]. The norm-based entanglement measures constructed from these two maps, 
called negativity and <&-negativity, respectively, lead to two sets of bounds on the entangle- 
ment of formation, the tangle, and the concurrence. We compare these bounds and identify 
the sets of 4 x N density operators for which the bounds from one constraint are better than 
the bounds from the other. In the process, we present a new derivation of the already known 
bound on the concurrence based on the negativity. We compute new bounds on the three 
measures of entanglement using both the constraints simultaneously. We demonstrate how 
such doubly constrained bounds can be constructed. We discuss extensions of our results to 
bipartite states of higher dimensions and with more than two constraints. 
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I. INTRODUCTION 

Characterizing quantum entanglement [JJ, [2] is an important open problem in quantum in- 
formation theory |3j - The nonclassical correlations associated with entanglement have been of 
immense interest since the very inception of quantum mechanics [3J, [5]. Quantum information 
science has identified entanglement as a potential resource. The ability of quantum computers to 
solve classically hard problems efficiently, the increased security of quantum cryptographic pro- 
tocols, the enhanced capacity of quantum channels — all these are attributed to entanglement [3] . 
The presence of entanglement has been related to quantum phase transitions and the behavior of 
condensed systems [6], [7], [8]. Entanglement has also allowed the understanding of techniques such 
as density-matrix-renormalization group in a new light [SJ. A significant part of recent research 
in theoretical quantum information science has centered around understanding and characterizing 
entanglement. In spite of this, entanglement remains a poorly understood feature of quantum 
systems. 

Although many tests have been devised which attempt to decide whether a general quantum 
state is separable or not, this problem is known to be NP-Hard [10 . Quantifying entanglement 
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involves devising functions acting on quantum states that, in some reasonable way, order entangled 
states according to the degree of nonclassical correlation possessed by them. Measures of entangle- 
ment can be broadly divided into two classes depending on whether an efficient way of computing 
them for arbitrary states exists or not. Tests for separability can also be classified in a similar 
fashion [2] . Computationally operational measures of entanglement are easy to calculate for any 
state, while there is no known procedure for efficiently calculating computationally nonoperational 
measures for an arbitrary state. From here on we abbreviate the descriptions computationally 
operational and computationally nonoperational to simply operational and nonoperational, respec- 
tively. Several physically significant measures of entanglement are of the nonoperational variety. 
This makes it important to place bounds on the values of such measures. In this paper, we in- 
vestigate the problem of placing lower bounds on nonoperational measures of entanglement for a 
quantum state assuming that we know the values of one or more operational measures for that 
state. 

The outline of this paper is as follows. In Sec. [TT] we start with examples of both operational 
and nonoperational measures of entanglement. We then discuss the general scheme of placing 
bounds on nonoperational measures using operational ones as constraints. In Sec. |III| we start with 
the separability criterion due to Breuer [11] and then show that a new, operational entanglement 
measure, called the ^-negativity, can be extracted from it. In Sec. |IV| we use the ^-negativity to 
bound three nonoperational measures of entanglement for 4x# systems, namely, the entanglement 
of formation, the tangle, and the concurrence. We compare our results to the bounds based on 
another operational measure, the negativity. In the process, we present a different way of deriving 
the results in [12] . In Sec. |V]we obtain bounds on the three nonoperational measures using both the 
negativity and ^-negativity simultaneously as constraints. We also discuss how our new bounds 
relate to previously known bounds in this section. Our conclusions and future prospects are 
summarized in Sec. IVII 

II. GENERAL CONSIDERATIONS 
A. Operational and nonoperational measures of entanglement 

A commonly used measure of entanglement for a pure-state |\&) of two systems A and B is the 
entropy of the reduced density operator pa (or ps), 

S(p A ) = -Tr(p A logp A ) = S(p B ) = -Tr(p B logp B ). (2.1) 

We write this entropy either as a function h(^f) of the state 1^) or as a function H{^i) of the vector 
of Schmidt coefficients of \^f). It is a physically motivated quantity, in that it gives the rate at which 
copies of a pure state can be converted, by using only local operations and classical communication 
(LOCC), into copies of maximally entangled states and vice versa [E]. This measure can be elevated 
so that it applies to bipartite mixed states also by taking the so-called convex-roof extension of 
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Eq. (2.1 ). This extended quantity is the entanglement of formation (EOF), and it is defined as 

h(p) = min < > p-hity 



p=Y / pj\*Mn\- ( 2 - 2 ) 



The EOF provides an upper bound on the rate at which maximally entangled states can be distilled 
from p and a lower bound on the rate at which maximally entangled states must be supplied to 
create copies of p |14j . Exact expressions for the EOF of several classes of states are known. One of 
the earliest, and simplest, was for an arbitrary state of two qubits [15]. The EOF in that case, was 
presented in terms of the concurrence, a subsidiary quantity. The concurrence itself has since been 
identified as an entanglement monotone and extended to higher-dimensional systems |16j,[17 . 

The EOF and the concurrence are examples of a more general framework of defining entan- 
glement measures. Suppose we have an entanglement measure g defined only on pure states |\&), 
which is a concave function G of Schmidt coefficients /x of the marginal density operator of \^f). 
That is, suppose g has the form g(^) = G(/x) on pure states. This can be extended to a measure 
on mixed states via the convex-roof extension, 

g{p) = min j P = £>;I%>(*;|}- (2-3) 

It has been proven [TH] that any g{p) constructed in this way is, on average, nonincreasing under 
LOCCs. An entanglement measure with this property is known as an entanglement monotone. 
Besides the EOF and concurrence, other examples of entanglement monotones include the tangle, 
relative entropy, entanglement of distillation, etc. Each has its use in particular physical contexts. 
All the entanglement measures just mentioned have one feature in common: they are nonopera- 
tional. The bottleneck in evaluating most of these measures for mixed states is the minimization 
over all pure-state decompositions. As a consequence, placing lower bounds on these measures of 
entanglement for arbitrary states becomes important. 

An alternate approach to detecting and quantifying entanglement is based on the application 
of positive (but not completely positive) maps on density operators [19 ,[20 , [21], [22], [23], [23], 
|25j . [26], [27], [28]. In particular, a quantum state is separable if and only if it remains positive 
semidefinite under the action of any positive map. Given a positive map, we can construct an 
entanglement measure based on the spectrum of the density operators under the action of the map 
|29j . [30 . Such measures are typically much easier to calculate for general quantum states than 
the ones discussed earlier because they do not involve the convex-roof construction. Measures of 
entanglement based on positive maps are therefore operational in nature. The negativity |31j . |30j 
is an example of an entanglement measure of this sort, derived from the transpose map [32], [24 . 

We can use the operational measures of entanglement as constraints to obtain bounds on nonop- 



erational, convex-roof-extended ones. The complexity of the minimization in Eq. (2.3) is reduced 
by solving it over a constrained set, instead of over all pure-state decompositions. This was done 
in [33], [T2] for the EOF and the concurrence by minimizing over states with a given value of 
negativity. We turn now to describing the general procedure for constructing bounds based on the 
use of one or more operational entanglement measures as constraints. 
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B. Multiply-constrained bounds on nonoperational measures of entanglement 

Let fx, ■ ■ ■ , fx be K operational measures used to characterize the entanglement in a bipartite 
system. Assume that they have values n = (m, . . . , nx) for a state p. Their action on pure states 
can be expressed as functions of the Schmidt coefficients, i.e., 

i = !,-■■ ,K. (2.4) 

We are interested in a lower bound on the value of an independent, nonoperational measure g, 
which is a monotone defined on mixed states via the convex-roof construction. Let us assume that 
for the state p, the optimal pure-state decomposition with respect to g is p = X^Pil^X^I- Then 

g(p)=J2Pj9^ j )=J2^ G (^)- ( 2 - 5 ) 

3 3 

Now define the function 

G(mi, . . .,m K ) = G(m) = min (G(/x)|Fi(/Lt) = mi, . . .,F K (n) = m K }. (2.6) 

Let <?(m) = co[G(m)] be the convex hull of G(m), i.e., the largest convex function of K variables 



(mi, . . . , mx) that is bounded from above by G(m). Using Eq. (2.6) and the convexity of Q, we 
can write 

g{p) > X>e(n') > sfepjiA (2.7) 

3 3 

If Q is a monotonically nondecreasing function of all its arguments and if the operational measures 
Fi are convex functions so that ^2jPjnj > m, we obtain 

g(p) > Q(n). (2.8) 



If the conditions for the validity of the inequality (2.8) are met, then we obtain a lower bound on 
g{p) by knowing the operational measures n for p. 



Regrettably, the first assumption leading to inequality (2.8) is not always valid: the function 
G(n) is not guaranteed to be monotonic. If it is not, then we have to impose monotonicity by 
introducing a new monotonically nondecreasing function G+(n), constructed from G(n). In the 
examples we consider in Sec.|vJ G(n) turns out to be monotonic, so we do not have to construct the 
new function Grt(n). For the sake of completeness, the general construction of Gj(n) is presented 
in Appendix [Aj 

We can now redefine £7(n) as the convex hull of Grf(n), rather than simply the convex hull of 
G(n). It is not immediately obvious that the convex hull of a monotonically nondecreasing function 
is also monotonically nondecreasing. The proof that this is so is given in Appendix [B] 

The only requirement on the operational entanglement measures Fi for using them as constraints 
is that they are convex functions on the set of states. Furthermore, even if we do not have 
the functions Fi themselves, but have instead functions Fi that bound Fi from above for pure 
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states, then the functions F{ can be used as constraints in the definition (2.6) of G, in place of 
the functions F{. The arguments leading to inequality (2.8) go through exactly as before, i.e., 



dip) — Y^jPjGi™ 3 ) — Q [Ylj Pj™ 1 j — @ylLjPj n ' > j — ^( n )) the only difference being that there is 
an additional step, the second-to-last one, where we use Fi(/j.) = hi > rii = -F, to conclude that 
Ylj Pj'H — Ylj Pj n j- The danger in using upper bounds instead of the actual values of the functions 
Fi is that the final bound on g(p) might turn out to be less useful or even meaningless. In the 
example we consider in Sec. |VJ however, one of the constraints we use is an upper bound on an 
operational entanglement measure, rather than the measure itself, yet the bound we get turns out 
to be stronger than previous bounds. 

Since our bound is intended for arbitrary states, there is one more subtlety to address, and that 
is the domain of the functions G, G, and Q. The operational measures n map the state p to a 
point in a i<C-dimensional hypercube in the space of the K independent constraints n^. Pure states 
correspond to a simply connected subset in this hypercube, which we call the pure-state region. 
The pure-state region is the domain of the functions G and G. This domain is not always convex, 
and so Q{vl) is defined on the convex hull of the pure-state region, which is generally bigger than 
the pure-state region, though only a subset of the full hypercube available to a general state. 

Finally, we have to extend <5(n) to the entire hypercube of states. Note that for inequalities (2.7) 
and (2.8) to hold, £/(n) must be a monotonically nondecreasing function in the entire hypercube 
while it has to be convex only on the convex hull of the pure-state region. So, in extending £7(n) 
outside the hull, we only have to take into account the monotonicity requirement (2.8 ). To construct 
such an extension of £7(n), start from a point on the boundary of the hull and begin traversing 
out along decreasing directions parallel to the axes of the hypercube. Outside the hull, and till 
reaching the boundaries of the hypercube, the extension is defined as the constant function with 
value equal to that at the point on the boundary of the hull. To generate the complete extension, 
this simple procedure is repeated for every point on all the boundaries of the hull. This procedure 
is also demonstrated in detail in Sec. |V]for the example we consider. 

In this paper, we carry out the general program just described with two particular constraints 
{K = 2). One of them is the negativity [30]. For the second, we develop a new entanglement 
measure, called ^-negativity, based on a recently presented separability criterion [11] (see |34j for 
another measure based on the same criterion) . Like the negativity, it is easily computable for any p 
and there are no convex-roof constructions involved in the computation. The ^-negativity, unlike 
the negativity, is not a simple function of the Schmidt coefficients for pure states. We find a simple 
function of the Schmidt coefficients that is an upper bound on the ^-negativity and, as described 
above, we use this function instead as the constraint to simplify our computations. We use both 
the (upper bound on) ^-negativity and the negativity simultaneously as constraints to place new 
bounds on the EOF, tangle, and concurrence of 4 x JV systems. Ours is the first instance of a 
doubly-constrained bound on entanglement measures for a family of states. It puts bounds that 
are tighter than those obtained in [33], [p2]. Multiply constrained bounds based on entanglement 
witnesses that can be applied to individual quantum states have been obtained using a different 
approach in [35 , [36 . 
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Although all of the results in this paper are obtained using the negativity and ^-negativity, 
a third constraint based on the realignment criterion (26], [27], [28] can be added to improve the 
bounds for certain classes of states. On pure states, the negativity and the realignment criterion 
lead to the same constraint. This means that in deriving both the singly and doubly constrained 
bounds we could have modified the negativity to take advantage of this, as was done in [33], [51] . 
Furthermore, the addition of the realignment criterion adds very little complexity to the procedure 
described below. 

Before concluding this section, we review the notation used in this paper. We use lower case 
Latin letters, say g, to denote entanglement measures. The corresponding upper case character, 
G, denotes the same entanglement measure defined on pure states, expressed as a function of the 
Schmidt coefficients. The same letter with a tilde on top, G, stands for the minimum of G subject 
to constraints. Calligraphic letters like Q denote the bound on g obtained by taking the convex hull 
of G. If we have to impose monotonicity on G as an intermediate step, we define a new function Gf. 

III. $-MAP 

Recently, a new separability criterion has been proposed based on a positive nondecomposable 
map [11. It is a combination of the Peres criterion and the reduction criterion |37| . |38] for 
detecting entangled states. In this section we construct a new entanglement measure from this 
map and calculate it for pure states. 

A. Separability criterion 

Let us consider a finite-dimensional Hilbert space C D . It can be regarded as the space of a 
spin-j particle with D = 2j + 1. A natural basis for this space is the "angular-momentum basis" 
| j, m), where m = — j, — j + 1, . . . , j — 1, j. The separability criterion to be presented involves the 
time-reversal operator i? whose action on an operator a acting on C D is given as 

$a = Va T V\ (3.1) 

where the superscript T stands for transposition in the angular-momentum basis and V is a unitary 
operator defined as 

(j,m\V\j,m') = (-iy~ m 5 m ^ m , (3.2) 

This map was initially introduced by Breuer to study the entanglement of 4 x 4 SU(2) invariant 
states; in that case, the $ map, together with the Peres criterion, was found to be a necessary and 
sufficient separability condition |39j . In even dimensions, an additional property holds: V T = — V, 
i.e., V is skew-symmetric in addition to being unitary. 

The condition for positivity under the partial time-reversal map (I ® $)p > is unitarily 
equivalent to the Peres PPT criterion (I ®T)p > 0. This means that partial time reversal can be 
used as an entanglement detection criterion. Breuer |11] defines a positive map 

$(p) = Tr(p)/ - p - Vp T V^ , (3.3) 
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which conjoins the time reversal map with the so-called reduction criterion [37]. The map $ then 
defines for any joint density operator pab a necessary condition for separability as 

(I ® $)p AB = Ti B ( PAB ) ®I B - p AB - (I A ® V)p T A B B {I A ® V ] ) > 0. (3.4) 

Any state that violates the above condition must be entangled. 

Consider the space can be regarded, without loss of generality, as the 

Hilbert space of two spin-j particles with j = {D — l)/2. The total spin of the system, denoted 
by J ranges over the values J = 0, 1, . . . , 2j = D — 1. Let Pj be the projector onto the (2J + 1)- 
dimensional spin- J manifold. It can then be shown that $ is a nondecomposable positive, but not 
completely positive map |34] in all even dimensions D greater than or equal to 4. The proof 
of positivity cannot be extended to odd dimensions as it exploits the skew-symmetric nature of the 
unitary operator V . In addition, the hermitian operator 

W = {I0^)P O (3.5) 

is an optimal entanglement witness (ID], in that the set of PPT states detected by W is not 
contained in the set detected by any other single witness. There, of course, exist families of PPT 
states that W fails to detect. The optimal nature of W provides motivation for contructing a 
measure of entanglement based on the $-map. 



B. Entanglement measures from maps 



Our endeavor here is to define a quantitative operational measure of entanglement based on the 
map. We call this quantity the ^-negativity, denote it by n$, and define it for a general mixed 
state as 



n<s>(p) 



D{D -I) 



(7® 
D-2 



1 



(3.6) 



where D = min(dim(H J 4), dim(HB)), and the trace norm of an operator is defined as ||0|| = 
For a separable state a, (I (g> $)<r has no negative eigenvalues, so ||(J (8) $)<r|| = 
Tr[(/(g) <J>)cr] = (D — 2)Tr(<r) = D — 2. Hence the ^-negativity is zero on separable states. This 
calculation also shows that $ is not trace preserving, and this motivates the factor of D — 2 in 



the denominator of Eq. (3.6). The $- negativity is a shifted and scaled version of the sum of the 



negative eigenvalues of a state under the action of the map (3.4). Since this sum can be expressed 
in terms of the trace norm of an operator, 1 1 (/ <8> &)pab 1 1 , it is a convex function of p as required in 
the general scheme described in Sec. [TTJ By defining the ^-negativity in terms of a map, we make 
sure that it is an operational measure that involves no convex-roof extensions. 

Similar measures of entanglement based on other positive, but not completely positive maps 
have previously been proposed and investigated [29 . The negativity [3D], which is based on the 
Peres partial transpose criterion, is defined as 

\\n T A\ 



n T (p) 



1 



(3.7) 
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where p is a joint density operator, Ta is the partial transposition with respect to system A. A 
positive value of nj- indicates an entangled state. The Peres negativity, in addition to being a 
measure of entanglement, is also an entanglement monotone, since it is nonincreasing on average 
under LOCC operations |30| . The ^-negativity is not an entanglement monotone, but it is a convex 
function of p. 

The ^-negativity is a new operational measure of entanglement for any quantum state. To use 
the ^-negativity as a constraint in bounding nonoperational measures we need expressions for n$ 
for pure states. We start from the Schmidt decomposition of any pure state, 

D 

\VAB) = y\y/m\ai,bi) (3.8) 



1=1 

for \^ab) 6 C D ® and D < N. The /i, are the Schmidt coefficients, satisfying pi > V i and 
Before we apply the map (/ (g) 3>) to this state we note that the matrix V appearing in the 



definition of the $-map has the form given in Eq. (3.2) only in the angular- momentum basis for 
system B, and the required transposition is also carried out in this basis. Relabeling the angular- 
momentum eigenvectors {\j,m)} as {|/)} with I = 1, . . . , D = 2j + 1, we transform I^ab) to the 
angular-momentum basis for subsystem B and obtain 

D 

Pab = \^ab){^ab\ = ^2 V^iVj {l\k}(bj\m}\ai,l)(aj,m\, (3.9) 

i,j,l,m=l 

and 

D D 

(I<g><S>)p AB = ^2tH\ai)(ai\<»^2\l)(l\ 
i=i i=i 

D 

~ XT ^PiPj{l\b i ){b j \m)\a i ,l)(a j ,m\ 

i,j,l,m=l 
D 

- ^MH{-^ l+m (b l \l)(m\b j )\a l ,D + l-m)(a j ,D + l-l\. (3.10) 

i,j,l,m=l 

The trace norm 11(1® /Oa_b 1 1 and the entanglement measure n$ defined using the trace norm are 
rather complicated functions of the Schmidt coefficients p% and the matrix elements {l\b%) of the 
unitary matrix that transforms between the Schmidt basis of subsystem B the angular momentum 



basis used in Eq. (3.2). Computing the numerical value of n$ for any state is relatively easy, but 
the analytic expression for the entanglement measure is quite unwieldy. 

All we really need to generate a constraint from the 3>-map, which can be used to place lower 
bounds on nonoperational measures of entanglement, is an upper bound on Such a bound is 
obtained by considering the special case in which the Schmidt basis of subsystem B is the same as 
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the basis of the angular-momentum eigenstates {\l}}, i.e., (l\bi 

D 

(I®$)p AB = y~]fJLj\ai)(a>i\ ® y ]\bj)(bj\ 



5u. We then have 



i=i 



D 

£ 

3=1 



D 

»J=1 



|aA)(aA'| + (— l)* + - 7 '|a i ){a i | ® |& i3 _ 7 - +1 )(&£)_ i+ i 



(3.11) 



For the first nontrivial case, D = 4, which we will be using extensively, explicit diagonalization of 
the above operator is possible. As shown in Appendix [C] it has six nonzero eigenvalues, one of 
which is negative. The trace norm can then be evaluated as the sum of the absolute values of the 
eigenvalues. Thus, 

1 1 (J® $)pab\\ = ® <5>) P AB?) = 2[1 + Vi^i + A»4)(/X2 + #0], (3.12) 

where we use the fact that the $-map is hermiticity preserving. Therefore, for all 4 x iV pure states 
that have the Schmidt basis for subsystem B the same as the angular-momentum basis, 



3\/ (Mi + Vi){H2 + M3) = n$. 



(3.13) 



In Eqs. (3.12) and (3.13) and in all our subsequent discussion of the function n$, the Schmidt 



coefficients are ordered from largest to smallest, i.e., p\ > fj>2 > P3 > M4- 

The function n$ is a simple function of the Schmidt coefficients for any 4 x N pure state and 
as shown numerically by the results displayed in Fig. [TJ the true $- negativity, n$, calculated with 
respect to a fixed angular-momentum basis, is always bounded from above by h^. 



Counts per bin 
1x10 

8x10 

6x10 

4x10 



2x10 




0.4 0.5 



no -n$ 



FIG. 1: A histogram of the difference n$, — n$ for five million 4x4 pure states picked randomly from the 
Haar measure. The bin size in the histogram is 0.001. The difference is always found to be positive. These 
results carry over to 4 x N states when N > 4 (N even), since the difference only depends on the Schmidt 
coefficients. We have tried to prove that the difference is nonnegative without success and thus rely on this 
numerical demonstration instead. 

In the rest of this paper we use n$ instead of n$ as the constraint for bounding nonoperational 
measures because of its simple algebraic form. When we refer to constraints based on the $- 
negativity we are referring to fixing the value of n$. Expressions for n<j> for pure states in higher 
dimensions are discussed in Appendix [Cj 
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In the next section, we will use n<j> to put lower bounds on the EOF, tangle, and concurrence for 
4 x N mixed states. We then compare our results to such bounds that have already been derived 
based on the Peres negativity tit- For D x N pure states, the Peres negativity is given by [31,[30j 




(3.14) 



IV. SINGLY CONSTRAINED BOUNDS 



A. Entanglement of formation 

A lower bound W(fi$) on the EOF, constrained by pure states having a certain value for n$, 
can be obtained using the steps described in Sec. [IT| All the subsequent results presented in this 
section and the next are for 4 x N states p, with N > 4. 

Firstly, we have to find 

£T(n$) = min |jJ(/x) 3^ (pi + M4XM2 + ^3) = , (4-1) 
and then its convex hull, 

H(p) = co[H(n*)], (4.2) 

provided if (n$) is a monotonically increasing function of n$. Defining pi+p^ = a and p 2 +P3 = 0, 
we can write the normalization and n$ constraints as 



ft 2 



a + (3 = I and ap=-£-, (4.3) 



which give 



l±Jl-4n|/9 l + Jl-4hl/9 
a = 1— and (3 = (4.4) 



Minimizing 

H(fjL) = -p 1 logpi - /i 4 log/i 4 - /U 2 log/i2 - /U 3 log/i 3 = H 2 (a) + aH 2 (pi/a) + (3H 2 (p,2/P), (4.5) 
where H 2 {-) is the binary entropy function, is trivial, because we simply make the last two terms 



zero by choosing p,\ = a and p 2 = j3 [and choosing the upper sign in Eq. (4.4 ) so as to be consistent 



with the assumed ordering of the Schmidt coefficients] . Then the minimum entropy is 

H{n*) =H 2 (a). (4.6) 

That H(h§) is a convex, monotonically increasing function of can be shown by considering 
its first and second derivatives. Its convex roof is the function itself, i.e., 

W(n«) = co[H(n*)] = H(m), (4.7) 



11 




and the bound can thus be extended to mixed states, giving 

l + ^/l-4n|/9' 

h{ P ) > H 2 • 



(4.8) 



with n<j> being the ^-negativity of p. 

The first step in bounding the EOF with only a single constraint on the negativity is to determine 
the function 

2 



Hlnj-) = min < 



H(jm) 



1 



riT 



(4.9) 



This was solved in [H],[33] for 2 or 3 Schmidt coefficients and recently shown to be valid for any 
number of Schmidt coefficients |42| . In particular, for four Schmidt coefficients, the case of interest 
here, we obtain 



H(n T ) =# 2 ( 7 ) + (l- 7 )log 2 3, 



(4.10) 



with 



7 



V2n T + 1 + ^3(3 - 2n T ) 
16 



(4.11) 



Unlike H(n§) , H[ht) is not convex over the entire range of nr. It is, however, a monotonically 
increasing function of ny. The actual bound on the EOF is thus the convex-roof extension of this 
function, co[.ff (nr)], which is given as [33] 



h{p) > H(n T ) = co [H(n T )] 



^2(7) + (1-7) log 2 3, n T e [0,1], 
(n T - §)log 2 3 + 2, nr €[!,§]. 



(4.12) 
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FIG. 3: In Region 1, the singly-constrained n$ bound is better than the singly-constrained ut bound. In 
Region 2, the opposite is true. 

Both the singly constrained bounds are plotted in Fig. |2j It might seem that the bound based 
on the ra$ constraint is always poorer than that in Eq. ( |4.12[ ), but this is not the case. There is a 



region in the h^-riT plane where the bound of Eq. (4.8) is better than that of Eq. (4.12). This is 
depicted in Fig. [3] 



B. Tangle and concurrence 

The procedure in the previous section can be undertaken for the tangle t{p) and the concurrence 
c(p) PS]) [HI- To place bounds on the tangle, we start by finding 



T(fi$) = min < 2 (l — \fi\ 



%V (Mi + M4)(A*2 + A*3) = ™* 



(4.13) 



which gives a bound for pure states. Then, just as for the EOF, the bound on the tangle for mixed 
states is given by the convex hull of T(re$), 

t{p) >T(n$) = co [?(n*)], (4.14) 

provided T(n$) is a monotonically nondecreasing function of h§. 



Using the normalization and n$ constraints of Eq. (4.4), we have 

/ 4 \ 



2(1-|H 2 )=2(1-J>?] =4^W 

i<j 



i=l 



n 3> 

i I — + M1M4 + 



(4.15) 



Just as for the EOF, the minimization is trivial, the minimum occurring for the upper sign in 

= (a* 2 = 

f = (4.16) 



Eq. (4.4), with \x± = = a) and fj,^ = (/i2 = /?), thus giving 

4.o 
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Since this is both monotonically increasing and convex in n<j>, the same bound holds for mixed 
states, but in terms of the actual negativity n$, i.e., 

4 



t{p) >T(n$) 



r> 2 



(4.17) 



The lower bound on the tangle, subject to a constraint on the negativity, is found by starting 
from 

( ' , , 777 



T{nr) 



mm < 



2(1-|H 2 ) 



1 



tit 



(4.18) 



This is a relatively involved minimization, but it is exactly the same as the minimization problem 
that arises in evaluating a bound on the tangle for isotropic states, so we can adapt the result 
of [T7] to give 

1 



T{n T ) 
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+ 4n^ + J 3 (3 + 4n T - 4n|) (2n T - 3) . 



(4.19) 



This quantity is monotonically increasing, but is not convex over the complete range of nx- The 
convex hull T(jit) = co [T(n^)] is required to extend the bound to mixed states. Again using the 
results of [17] , we obtain 



t( P ) > T(n T ) 



\(9 + An 2 T + yj3 (3 + 4n T - 4n 2 ) (2n T - 3) J , n T E [0, 1] , 

n T E [1,1]. 



i 

2' 



(4.20) 



We can derive from Eq. (4.17) an expression for the lower bound on the concurrence of 4 x N 



states with a given value of n$: 



z{p) > C(n$) = C(n^,) = \JT{n^) = -n$. 



(4.21) 



An expression for the minimum of the concurrence, subject to the negativity constraint, can be 



obtained from Eq. (4.19). The resulting function is everywhere concave, and thus its convex hull 
is a straight line joining the end points. This line is 



c(p) > C(n T ) = y|^T- 



(4.22) 



The bounds on both the tangle and the concurrence are plotted in Fig [4] As was true for the 
EOF, the n$ bound is better than the ut bound in some parts of the n^-riT plane. This is shown 
in Fig. [5] 

Recently, a lower bound on the concurrence has been derived based on the negativity con- 
straint [12., using techniques different from those employed here. That lower bound is exactly the 



one in Eq. (4.22). We have thus provided an independent derivation of the bound presented in [12 . 
In addition, we can use the procedure from [12] to derive a lower bound on the tangle based on 
the n$ constraint. Then we obtain 



T(n$J 



(4.23) 



which for general mixed states, leads exactly to the bound in Eq. (4.17). 
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FIG. 4: The plot on the left shows the bounds on the tangle and the concurrence based on the n$ constraint. 
The solid line is the bound on the tangle and the dashed line is the bound on the concurrence. On the right 
is a plot of the analogous bounds based on the jit constraint. 





FIG. 5: Region 1 is where the n$ constraint is better than the tit constraint for bounding the tangle and 
concurrence. Region 2 is where the converse is true. 

V. DOUBLY CONSTRAINED BOUNDS 



In this section we place new lower bounds on the EOF, tangle, and concurrence for 4 x N density 
operators by using ut and n$ simultaneously as constraints. 
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A. Pure states of 4 x N systems 



For a4xJV pure state, described by the Schmidt coefficients fn, i 
constraint equations, 



1 , . . . , 4, we have three 



3V (Ml + M4)(yW2 + A*3) = n<f>, 

Ml + M2 + M3 + M4 = 1, 



in addition to the inequality constraints 



1 > Mi > M2 > M3 > M4 > 0. 



(5.1) 



(5.2) 



Both n$ and take on values between and 3/2, so all 4 x iV states, pure or mixed, are mapped 
to a square of side 3/2 in the re^-n-r plane. Not all points in the square correspond to pure states. 



If we solve the three equations in (5.1 ) simultaneously and express m, \ii and /X3 in terms of n^, h§ 



and fi4 (see Appendix [D]) , we find that for some allowed values of n$ and n-p, there is no allowed 
value of (14 for which the other three Schmidt coefficients are real numbers between and 1 in even 



one of the solution branches of (5.1 ) 



To find the region occupied by pure states in the n$-«7 plane, let us use the pure-state expres- 



sions for riT and in Eq. (5.1 ) to find the largest and smallest values that tit can take on for a 
fixed value of n$. We proceed as in the minimization of H(fi) in Sec. IV Defining a = \x\ + /U4 and 
(3 = M2 + M3> the normalization and n$ constraints can be solved to give a and (3 as in Eq. (4.4). 
The negativity takes the form 



yj2n T + 1 = v^i + v 7 " - Mi + \/M2 + V/3- M2 • 



(5.3) 



It is trivial to see that the maximum of nx occurs when [i\ = ^4 = a/2 and Hi = ^ = f3/2. This 
maximum cannot be achieved, however, because we must respect the ordering /^i > /i2 > M3 — M4 
that we assumed in our definition of n<j>. We should always choose \i2 = M3 = /?/2, but the best 
we can then do with [i\ and ^4 is to choose ^4 = /3/2, ^1 = a — (3/2 when a > (3 [upper sign 
in Eq. (4.4)] or fi\ = (3/2, ^4 = a — (3/2 when (3 > a [lower sign in Eq. (4.4)]. The requirement 
that nx < a implies that the latter case can only be used when n$ > y/2. In both cases, the the 
maximum value of nj< for fixed n<j> has the form 



tit 



a - (3/2 + 37/372 -1 



(5.4) 



It turns out that the upper sign in Eq. (4.4) always gives a larger value for rvp- Using the upper 



sign, we find that the maximum of ut for fixed values of n$ is given by 



riT 




+ 2W1 



(5.5) 
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The minimum value of nx occurs on the boundary of allowed Schmidt coefficients, i.e., when 



[i\ = a and \i<i = (3, with the upper sign in Eq. (4.4). Thus the minimum value of tit for a fixed 
value of n$ is given by 



tit 



1 



1 



(5.6) 



From Eqs. (5.4) and (5.6) we find that the pure states of a 4 x N system lie in the region shown 



in Fig[6j Notice that for this case of two constraints, the pure-state region is not convex. 




0.2 0.4 0.6 0. 



FIG. 6: The pure-state region in the n$-nr plane for 4 x N systems. 



B. Entanglement of formation 



The EOF for pure bipartite states is a concave function of the marginal density operator obtained 
by tracing over one of the subsystems. This means that it is a concave function of the Schmidt 
coefficients fi. Searching for a minimum is not the most natural thing one can do with a concave 
function, yet this is what we are instructed to do by the procedure for bounding the EOF outlined 



in Sec. II B Starting from the EOF H([i) for pure bipartite 4 x N states, our objective is to find a 



convex, monotonic function H(n) as outlined in the Sec. [TTJ This function will be our lower bound 
on the EOF for all states. 

The first step is to find the function 

X.7 v^i' 



H(n) = H(h^,nr) = min < 



3 \/(mi + A*4)(a*2 + M3) = n$, 



1 



riT 



(5.7) 



which is defined on the pure-state region. The method of Lagrange multipliers is not suitable for 



finding the minimum in Eq. (5.7) because the problem is over-constrained. The equations that we 
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obtain using Lagrange multipliers have a consistent solution only if h$ and ut are related as in 



Eq. (5.5 1 and therefore lie on the upper boundary of the pure-state region. This does not mean 
that there is no minimum for H(fi), but rather that the minimum lies on a boundary of the allowed 
values of fx. 

We already know H(h§, rvA on the boundaries of the pure-state region. The boundary with 
three of the Schmidt coefficients being zero is the origin in the n$-nr plane where H(n§, tit) = 
H{ijl) = 0. The boundary with two of the Schmidt coefficients zero lies on the line = h§/3. To 
find the value of H(n§, tit) along this boundary, note that the minimum of the EOF subject to 
just the n$ constraint occurs for fjb$ = (a, 1 — a, 0,0), where a is given in Eq. (4.4). Substituting 
into tit we get ny = yj a(l — a) = n<j>/3. This means that along the line tit = n$/3, the ny 
constraint is automatically satisfied if the n$ constraint is satisfied. Thus along the lower boundary 
of the pure-state region, we have H{n^ : tit) = H(n$,). Similarly, along the upper boundary of 
the pure-state region, the n<j> constraint comes for free. This is because the minimum of the EOF 
subject to the tit constraint occurs when the Schmidt coefficients are given by fi T = (7, 7', 7', 7') 
with 7 given by Eq. (4.11) and 7' = (1 — 7)/3. The doubly-constrained problem reduces to the 
singly-constrained problem when n<j> = y/2(2^ + 1)(1 — 7). Relabelling 7 as a — (3/2 and 7' as f3/2 
we see that the n$ constraint is automatically satisfied along the upper boundary of the pure-state 
region if the tit constraint is satisfied. Hence along the upper boundary of the pure-state region, 
we have H(n§, tit) = H(nr). 

These considerations mean that for the entanglement of formation, the monotone boundaries 
that we define in Appendix [A] coincide with the boundaries of the pure-state region, making it 
unnecessary to construct the monotonically nondecreasing function (n$, tit) , since H{n^,riT) 
is itself monotonically nondecreasing. 

The minimum of H(fj.) in the remaining part of the pure-state region can be found using 
the straightforward numerical procedure described below. We start from the two distinct sets of 
solutions and of the three constraint equations (see Appendix [d|. We go to the boundary 
where one of the Schmidt coefficients is zero by setting ^4 = in the solutions. Now compute 
H(fj,W) and if (/x^ 2 **) corresponding to the two solutions in the regions in the ti^-tit plane where 
each of the solutions is valid. The solutions are not valid in the whole pure-state region because 
the three Schmidt coefficients have to be real, nonnegative numbers less than one. All points in the 
pure-state region cannot be covered if we set ^4 = 0. This is easily seen by noticing that the point 
n$ = riT = 3/2 corresponds to the fully entangled 4 x N state and for this state all four Schmidt 
coefficients have the value 1/4. The fully entangled state and other states close to it cannot be 
reached using the procedure described above if we stay on the boundary defined by /U4 = 0. So we 
start increasing the value of [i^ in small steps until it reaches 1/4. The parts of the 2-constraint 
region that are covered by different choices of \i± are shown in Fig. [7j 

This numerical procedure gives us ranges of values of /U4 over which H(n {1) ) and/or #(/i (2) ) 
can be calculated at each point in the pure-state region. For the value of H(n) at each point, we 
pick the minimum over the allowed range of values for at that point. 

The function H(n) in the pure-state region is shown in Fig. [s] It is, as required, a monotonically 
increasing function of both n<j> and tit- Along the upper boundary of the pure-state region, the 
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FIG. 7: The part of the 2-constraint region in which a value for H(n) can be computed is shown for four 
values of LI4 = 0, 0.02, 0.1, and 0.25. The two lines are the boundaries of the pure-state region. 



c 1 





FIG. 8: (Color online) Plots of H(n), the minimum of the entropy of formation, H(fj,), in the pure-state 
region. On the left side is a 3-dimensional plot of H(n) and on the right is a contour plot of the same 
function. 



numerically computed value of H(n) matches the value of H(nr) from Eq. (4.9 1 . In addition to 
this, from the contour plot of H(n) in Fig. ^ we see that along the upper boundary, the function 
has zero slope along the n$ direction. 

The function H(n) is not convex, which can be seen by computing the Hessian at every point 
in the pure-state region. If the function were convex, both eigenvalues of the Hessian would be 
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positive at all points. It turns out that one of the eigenvalues of the Hessian is negative in a region 
in the upper right corner of the h^-riT plane, close to the maximally entangled state. 
Since H(n) is not convex, we have to compute its convex hull, 



W(n) 



CO 



H(n) 



(5.8) 



to obtain the bound on the EOF in the pure-state region. The convex hull of H(n) can be computed 
numerically, and it turns out that the difference between H.(n) and H(n) is quite small (~ 10~ 3 ), 
the two differing differ only in a small region in the upper right corner of the pure-state region. As 
shown in Appendix [B] taking the convex hull preserves monotonicity. 





FIG. 9: (Color online) The doubly-constrained bound TC(n) on the EOF of all 4 x N states. On the right 
side is a contour plot of the same function. 



To obtain a bound on the EOF of all 4 x N states, we have to extend H(n) out of the pure-state 
region to the rest of the n$-nx plane. The extension has to respect the monotonicity of H(n) so 



that the string of inequalities Eq. (2.8) holds. This can be achieved by extending Tl(n) using 



surfaces that match the function at the lower and upper boundaries of the pure-state region. To 
preserve monotonicity, the surface added on to the region below the lower boundary of the set of 
pure states has zero slope along the ut direction, and the surface added on to the region above 
the upper boundary of the set of pure states has zero slope along the n<j> direction. The resulting 
doubly-constrained bound TL(n) on the EOF is shown in Fig. [9} We see from the figure that the 
extension to the whole n$-n^ plane produces a smooth and seamless surface. 

One final point worth mentioning involves the use of our bound for general mixed states. To 
do so, one must calculate n$ for the mixed state, and this calculation depends on the choice of 
an angular-momentum basis for system B in order to define the $-map. The bound itself thus 
depends on this choice of basis, and the best bound would generally be found for the basis choice 
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that gives the largest value of n$. For pure states, for example, the results in Fig. [I] show that the 
best choice of basis is the Schmidt basis for system B. 

The isotropic states, which lie along the diagonal in the h^-nx plane, are special in that they 



saturate the singly-constrained bound TL{tit) from Eq. (4.12 ). These states thus furnish a good con- 



sistency test of our doubly-constrained bound because our bound must match the singly-constrained 
bound when applied to isotropic states. A comparison of the two bounds for isotropic states is 
given in Fig. [T0| 




FIG. 10: The thick black line is the doubly-constrained bound on the EOF for isotropic states. The dashed 



white line, lying on top of the black line, is the singly-constrained bound H{iit) from Eq. (4.12) 



We can make a second comparison between the singly and doubly constrained bounds using 



Fig. 10 From the way we constructed 7i(n), we know that its value on the diagonal in the h^-riT 



plane is the same as its value on the upper boundary of the pure-state region. We also know that 
the upper boundary is where the singly-constrained bound and the doubly-constrained bound are 
the same. From Fig. 10 we see that the convex hull Ti^nx) of the function H(jit) of one variable 



matches the convex hull TC(n) of the function H(n) of two variables on the upper pure-state 
boundary. These consistency checks give us increased confidence in the accuracy of our results. 

C. Tangle and concurrence 



Doubly-constrained bounds can be placed on the tangle and the concurrence of 4 x N states by 
extending the procedure used for the EOF. For the tangle, we start by finding the function, 



T(n) = min < 



2(1- M 



3\/(A*i + M4X/U2 + M3) 



1 



tit 



(5.9) 



in the pure-state region. For the concurrence, we want the function C(n) = yT(n), since for pure 
states the concurrence is the square root of the tangle. 
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1.5 




FIG. 11: (Color online) The doubly-constrained bound T(n) on the tangle of 4 x N states. On the right 
side is a contour plot of the same function. 



1.5 




FIG. 12: (Color online) The doubly-constrained bound C(n) on the concurrence of 4 x N states. On the 
right side is a contour plot of the same function. 

For all three of the entanglement monotones, EOF, tangle and concurrence, the monotone 
boundaries we define in Appendix |A| coincide with the boundaries of the pure-state region. This 
is because the singly-constrained bounds for all three measures correspond to the same sets of 
Schmidt coefficients, fi T = (7, 7', 7', 7') and jLt$ = (a, 1 — a, 0, 0), and we have already seen for the 
EOF that these Schmidt coefficients define the boundaries of the pure-state region. This makes 
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it unnecessary for these entanglement monotones to construct the monotonically nondecreasing 
function discussed in Appendix [A} In general, for two different measures of entanglement and two 
constraints, the singly-constrained bounds for the two measures need not correspond to the same 
Schmidt coefficients. 

Once we have T(n), the convex hull of this function extended to the whole h^-riT plane is the 
doubly-constrained bound on the tangle, T(n). A three-dimensional plot and a contour plot of 



T(n) are shown in Fig. 11 



The bound on the concurrence is the convex hull of the surface obtained from C(n). The 



resulting bound on the concurrence, C(n), is shown in Fig. 12 



VI. CONCLUSION 



We focused on two aspects of the problem of quantifying entanglement in this paper. The first 
was a comparison between the bounds on different measures of entanglement obtained by using 
ut and h§ independently as constraints. The second was the construction of doubly-constrained 
bounds on the three measures of entanglement that we considered. 

Starting from the <&-map [11], we found that we can define an entanglement measure, which 
we call the ^-negativity. The ^-negativity of arbitrary quantum states can be calculated in a 
straightforward manner, just like their negativity. We also found that we can obtain a much 
simpler function n$ of the Schmidt coefficients of pure states that is an upper bound on their 
^-negativity. Previous work |33|.|12j has shown that the negativity can be used as a constraint 
to place bounds on the EOF, the tangle, and the concurrence of bipartite states. We obtained a 
different set of bounds on these three measures of entanglement for 4 x N mixed states by using 
n$ instead as the constraint. The scheme for placing lower bounds on nonoperational measures of 
entanglement is general enough to allow us to use n$ instead of n<j> as the constraint. We were 
then able to compare the two sets of bounds on the measures of entanglement coming from using 
either one of the two operational entanglement measures as a single constraint. 

We found that the h^-nr plane for pure states can be divided into two regions depending on 
which constraint led to the better bound on a given measure of entanglement. This prompted us 
to consider whether we can construct a single, composite bound for each measure of entanglement, 
applicable to the entire h^-riT plane, by using both constraints simultaneously. It turned out that 
for 4 x N systems this is a tractable problem, and we obtained doubly-constrained lower bounds 
for the first time for the EOF, the tangle, and the concurrence. We showed how the bounds on the 
different measures of entanglement obtained for pure states can be extended to include all states. 
We found that the requirement of monotonicity on the bound defined on pure states dictates how 
to extend the bound to all states. 
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APPENDIX A: GENERAL CONSTRUCTION OF G T (n). 



In this Appendix, we describe the general procedure for constructing the monotonically nonde- 
creasing function Gf(n), which replaces G(n) when the latter function is not itself monotonically 
nondecreasing. 



As mentioned in Sec II B pure states of the system correspond to a simply connected subset, 
called the pure-state region, in the state hypercube in ~K K ; the function G(n) is defined only on 
the pure-state region. Within the pure-state region, we can define K hypersurfaces Sk as those 
on which the fcth constraint equation, Fk(fi) = n k, is automatically satisfied if the remaining 
K — 1 constraint equations are satisfied. We denote the value of n/% on Sk by n* k (n') where n' = 
(ni, . . . , nfc_i, nfc + i, . . . , uk)', the function n* k (n') can be regarded as the defining equation for Sk- 
On the hypersurfaces Sk, G(n) is effectively defined by K — 1 constraints. We denote the value of 
G(n) on Sk by Gk(n'). The minimum of any function subject to K constraints is always greater than 
or equal to its value when subject to K— 1 of these constraints, so we have G(n', rik) > Gfc(n'), where 
we have let (n',rik) = n. The inequality is saturated when rik = n* k (n'). Now consider G(n',nk) as 
a function of n&. If we fix n' and increase nk, starting from its lowest value, then G(n',n/ C ) has to 
either decrease or remain constant until we cross the hypersurface Sk- For n k > n k (n'), G(n',nk) 
is a nondecreasing function of n&. We want G-f(n) to be a nondecreasing function for all rik, so we 
define it by 

G^{n,n k ) = < ~ k = l,...,K. (Al 

I G(n) n k > n* k (n'), 

The construction of G|(n) is not complete at this point. Within each (K — l)-dimensional hy- 
persurface, we will encounter {K — 2)-dimensional hypersurfaces where two of the constraints are 
automatically satisfied. Across each of these {K — 2)-dimensional hypersurfaces, we can update 
the value of G-f(n) just as described above. 

There can be at most K different {K — l)-constraint regions and the if-constraint region will, 
in general, be surrounded by these {K — l)-constraint regions. The {K — l)-constraint regions are 
surrounded, in turn, by {K — 2)-constraint regions and so on. This construction procedure evi- 
dently terminates after K — 1 steps. We call the hypersurfaces identified in this appendix monotone 
boundaries, because the nested structure of /c-constraint regions they define are the key to con- 
structing the monotonically nondecreasing function G|(n) from G(n). In the examples we consider 
in Sec. |Vj the monotone boundaries coincide with the boundaries of the pure-state region, so we 
do not have to construct the function Gj(n). 



APPENDIX B: THE CONVEX HULL AND MONOTONICITY 



Here we show that the convex hull of a monotonically nondecreasing function on R K is also 
monotonically nondecreasing. We first define a partial order on the set of vectors in M. by defining 
x > y to mean Xk > Uk for all k. Define a monotone to be a function / : T> \— > [0, 1] satisfying the 
following conditions: 
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1. The domain D is a bounded region contained in the positive orthant (including boundaries) 
of R K , 

2. G V and /(0) = 0, 

3. V x, y £5), if x > y, then /(x) > /(y) . (monotonicity) 

The function / can alternatively be viewed as a set of points in R A+1 given by the tuples 
(xi, ... ,xk, f (*■)). Viewed this way, we can define the convex hull C of / as a set to be the 
smallest convex set containing the set /. We can also define the function c : T>' ^ [0, 1], to be the 
convex hull of / as a function. Thus c is the largest convex function bounded from above by /; in 
this paper c is called the convex roof of /; Clearly, c is just the lower boundary of the set C along 
the direction of the (K + l)st coordinate in 

Before continuing to the main theorem, we state an important result known as Caratheodory's 
theorem [43 . This theorem uses the notion of a generalized simplex of dimension d, which is just 
the convex hull of a set of d + 1 affinely independent points. A triangle, or example, regardless of 
shape, is a generalized simplex of dimension 2. For convenience we refer to a generalized simplex 
as just a simplex. 

Theorem 1 (Caratheodory) Let f be any bounded set of points in IR^ +1 , and let C = co[/] be 
the convex hull of f (as a set). Then x S C if and only i/x can be written as a convex combination 
of K + 2 (not necessarily distinct) points in f. Furthermore, C is the union of all the simplices 
with dimension less than or equal to K + 1 whose vertices belong to f. 

From Caratheodory's theorem and the fact that the function c is the boundary of the set 
C, we know that c can be expressed as the union of many simplices (usually infinitely many) 
whose vertices belong to /. These simplices are necessarily of dimension at most K, since the 
dimension of c is K. We can speak meaningfully about directional derivatives on these simplices 
and on c because of the following beautiful fact: any convex function has well defined one-sided 
directional derivatives everywhere and, furthermore, is differentiable everywhere except possibly a 
set of measure zero |43j . 

Theorem 2 Let f : T> \— ► [0, 1] be a monotone, and let c : T>' i— > [0, 1] be the convex roof of the 
function f . Then c is also a monotone. 

Proof: The domain V of c in general contains the domain T> of /, but it will remain bounded and 
in the positive orthant of M. K and is furthermore always convex even if T> is not. Clearly £ P', 
since G T>. The fact that c(0) = can be seen by the fact that /(0) = is the global minimum 
for /, and the convex hull of a function will always contain the function's global minimum. This 
shows that c satisfies the first two criteria of a monotone. 

Now we prove the final criterion, the monotonicity of c. Consider the set of all possible simplices 
with dimension less than or equal to K with vertices lying in /. From Caratheodory's theorem, c 
is a union of some subset of these simplices. However, every simplex in this set has the property 
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of monotonicity over its domain of definition. This follows from the "multidirectional" version of 
the mean value theorem jH] , |45 , for which we now sketch the proof. Suppose we choose a simplex 
s£c. Along a given direction p, the smallest value of the directional derivative of / lying above s is 
a lower bound on the directional derivative of s. In particular, if p > 0, then from the assumption 
of monotonicity of /, we know that V p / > everywhere, and hence V p s > 0. This implies that 
each constituent simplex in c is indeed monotonic. To show that c is a monotone, we use the 
convexity of c to see that the directional derivative in some direction p > across two neighboring 
simplices s% and S2 cannot decrease. 



APPENDIX C: n<s> FOR D x N PURE STATES 

Our objective in this Appendix is to characterize the eigenvalues of the operator (I®&)pAB = O 



in Eq. (3.11 ) for the special case in which the Schmidt basis for subsystem B of the pure state pab 



is the same as the angular-momentum basis. Recall that in Eq. (3.11), the Schmidt coefficients 
Hi are ordered from largest to smallest. The operator O is a DN x DN operator, although it 
clearly has rank at most D 2 , so we can regard it as a D 2 x D 2 operator, having D 2 eigenvalues. 
Although O can be written in matrix form in the Schmidt basis, we refrain from doing so here, 
as the expression is unwieldy and not very illuminating. We can, however, by permuting the rows 
and columns of O, write it as 

O = O0TeR, (CI) 

where is a matrix of zeros, of size D x D. 

To describe T and R, we first make some definitions. An index is an integer between 1 and D. 
An ordered pair of indices (j, k) is said to be inadmissible if k = D — j + 1 or k = j. All other 
indices are said to be admissible . A set of indices is called admissible if the elements are pairwise 
admissible. A product of n distinct Schmidt coefficients Pj 1 Pj 2 ■ ■ ■ Pj n is said to be n-admissible if 
all of the indices are pairwise admissible, and if, in addition, j\ < j% < ... < j n - Finally, define S n 
as the sum over all n-admissible products. Then 



© WW (C2) 



admissible 

where each W(p ; q) is a 2 x 2 matrix of the form 

For each index, there are D — 2 other indices with which it can form an admissible pair. Hence, 
D indices form exactly D{D — 2)/2 distinct admissible pairs, and that is the number of possible 
W(p.<?)' s °f * ne given form. W( p ^ has eigenvalues and \i v + p q . Thus, T has D(D — 2)/2 zero 
eigenvalues and an equal number of eigenvalues \i v + fj, g , where (p, q) is an admissible pair. 
The matrix R has elements 

Rj/fc = -yJP-j^kiX - <^,fc)(l - 8j,D-k+l), (C4) 
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where j, k = 1, • • • , D. It is thus a D x D matrix. The characteristic polynomial of this matrix can 
be written as 

g{z) = z D ^ ^/ 2 + ^%(-l)*<S t+1 ^/ 2 -*- 1 - (Z?/2- l)(-l)^ 2 X^ , (C5) 

where 

D/2 

1= llbij+HD-j+i). (C6) 

It is evident that the matrix R has D/2 zero eigenvalues. The remaining eigenvalues are the zeroes 
of the function 

r D (z) = z D ' 2 + £ t{-lfS t+l z D ' w - [D/2 - 1)(-1) B / 2 X. (C7) 

i=0 

The Descartes rule of signs tells us that the above equation has no more than one negative root. 
In fact, if all the Schmidt coefficients are nonzero, there is exactly one negative eigenvalue, the 
negative root of ro{z). Otherwise, all the eigenvalues are nonnegative, and the pure state under 
consideration could be separable. 

Putting all this together, we conclude that the spectrum of O has 

1. D + D{D - 2)/2 + D/2 = D{D + l)/2 zero eigenvalues, 

2. D(D — 2)/2 positive eigenvalues of the form \x v + fj, q , where (p, q) is an admissible pair, and 
D/2 — 1 positive eigenvalues, which are the positive roots of Vd(z) = 0, 

3. One negative eigenvalue, the negative root of r£,{z) = 0. 

As an example, we present the case of D = 4. Then Eq. (C7| becomes r^(z) = z 2 — (/ii + 



Hi){H2 + fJs), which has zeroes ±\/ + ^a){^2 + M3)- 



For D = 6, the function (C7| is 

r 6 (z) = Z 3 - z(fJ,ifJL 2 + MlA*3 + M2/"3 + MlM4 + + + M3M5 + M4^5 

+ M2^6 + A*3^6 + M4M6 + + 2(^1 + /XeX/^ + A*5)(>3 + M4)- (C8) 

APPENDIX D: SOLUTIONS OF THE CONSTRAINT EQUATIONS 



The three constraint equations, 



3v(Mi + M4)(^2 + A*3) = w* 5 

A*i + ^2 + A*3 + M4 = 1, (Dl) 
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can be solved to express n\, fi2, and //3 in terms of ray, n$, and ^4. There are four sets of solutions, 
of which only two are distinct because the other two can be obtained by exchanging fi2 and /X3. 
The constraint equations are invariant under this exchange. The two distinct solutions are the 
following: 



and 



/4" 



4" 



M 2 » 



4 2) 



(2) 
M 3 



Here Go and 0i are given by 



l + A /l--n|-2 M4 



1 - 4 r? 2 
9 * 



1 - -r? 2 
9 * 



+ 2y / a -£i(/4 1) ) 

-2Jg q -Gi(/?) 



4 
9 ? 



1 + Wl- -n| + 2 W & -0i Mi' 



,( 2 ) 



1 + Wl- -hi 



2^0 -Gi^?)^ 



(D2) 



0o = 1 + 8(n T + ^4)v / M4(2n T + 1) - An T {n T + 4/i 4 ) - 3^ 4 (/Lt 4 + 2), 



(D3) 



(D4) 



£1(^1) = Y2 + (^) 7 (\Z2nT + 1 - VMi) + y (3 + 8n T - 8vW2nr + 1) + 5^4) 
! V ' ~ (l - 2v//i 4 (2n T + 1) + ^4) - n T {y/2n T + 1 - 3^mI) 



/i4 



(D5) 
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